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Summary. — The OPERA experiment in the underground Gran Sasso Laboratory 
(LNGS) has been designed to perform the first detection of neutrino osciUations in 
direct appearance mode in the muon neutrino to tau neutrino channel. The detector 
is hybrid, being made of an emulsion/lead target and of electronic detectors. It is 
placed in the CNGS neutrino beam 730 km away from the neutrino source. Runs 
with CNGS neutrinos were successfully carried out in 2008, 2009, and 2010. After 
a brief description of the beam and the experimental setup, we report on event 
analysis of a sample of events corresponding to 1.89 x 10^® p.o.t. in the CERN 
CNGS beam that yielded the observation of a first candidate Vt CC interaction. 
The topology and kinematics of this candidate event are described in detail. The 
background sources are explained and the significance of the candidate is assessed. 
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1. Introduction 

Two types of experimental methods can be used to detect neutrino oscillations: ob- 
serving the appearance of a neutrino flavour initially absent in the beam or measuring 
the disappearance rate of the initial flavour. In the latter case, one must know the flux 
of the beam precisely. In this type of experiment one explores whether less than the 
expected number of neutrinos of a produced flavour arrives at a detector or whether the 
spectral shape changes if observed at various distances from a source. Since the final 
state is not observed, disappearance experiments cannot tell into which flavor a neutrino 
has oscillated. An appearance experiment searches for possible new flavours of neutrino, 
which does not exist in the original beam, or for an enhancement of an existing neutrino 
flavour. The identification of the flavour relies on the detection of the corresponding 
lepton produced in its charged current (CC) interactions: viN — >■ l~X with I = e,fi,T 
and where X denotes the hadronic final state. 
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In the past two decades, several experiments carried out with atmospheric and accel- 
erator neutrinos, as well as with solar and reactor neutrinos, have established the picture 
of a three-neutrino oscillation scenario with two large mixing angles. Atmospheric sector 
flavor conversion was first established by the Super-Kamiokande [1] and MACRO [2] ex- 
periments and then confirmed by the K2K [3] and MINOS [4] longbaseline experiments. 
The CHOOZ [5] and Palo Verde [6] reactor experiments excluded indirectly the i^f^ — > i/g 
channel as the dominant process in the atmospheric sector. However, the direct obser- 
vation of flavour transition through the detection of the corresponding lepton has never 
been observed. Appearance of Ur will prove unambiguously that i/^- oscillation is 

the dominant transition channel at the atmospheric scale. 

The OPERA experiment [7] has been designed to directly observe the appearance of 
in a pure beam on an event by event basis. The signature is given by the decay 
topology and kinematics of the short lived r~ leptons produced in the interaction of 
Vt-N — > X and decaying to one prong (/i, e or hadron) or three prongs, which arc [8]: 

T- irvt,vr with BR = 17.36 ± 0.05% 
T- e-VeVr With BR = 17.85 ± 0.05% 

T- h- {m:^)Dr with BR ^ 49.52 ± 0.07% 
T~ -> 2h-h+{mT°)Dr with BR = 15.19 ± 0.08%. 

2. The Neutrino Beam 

The CNGS ly^ beam produced by the CERN-SPS is directed towards the OPERA 
detector, located in the Gran Sasso underground laboratory (LNGS) [9] in Italy, 730 
km away from the neutrino source at CERN. In order to study z^^ Vr oscillations in 
appearance mode as indicated in the atmospheric neutrino sector, the CERN Neutrinos 
to GranSasso (CNGS) neutrino beam [10] was designed and optimized by maximizing 
the number of Vt-GG interactions at the LNGS. 

The average beam energy is 17 GeV ^ well above tau production energy treshold. 
The Vfj, contamination is 4% in flux, 2.1% in terms of interactions. The Ue and J>e 
contaminations are lower than 1%, while the number of prompt Vr from Ds decay is 
neghgible. The average L/Ei, ratio is 43 km/GeV, suitable for oscillation studies at 
atmospheric Am^. Due to the Earth's curvature neutrinos from CERN enter the LNGS 
halls with an angle of about 3° with respect to the horizontal plane. 

With a nominal CNGS beam intensity of 4.5 x 10^^ protons on target {p.o.t.) per year, 
and assuming Ani^^ = 2.5 x 10~^eV'^ and full mixing, about 10 i/t events are expected 
to be observed in OPERA in 5 years of data taking, with selection criteria reducing the 
background to 0.75 events. 

The goal is to accumulate a statistics of neutrino interactions corrcspomding to 22.5 x 
10^9 p.o.t. in 5 years. The 2008, 2009 and 2010 runs achieved a total intensity of 1.78 x 
10""^^, 3.52 X 10^® and 4.04 x 10^^ p.o.t. respectively. Within these three years, neutrinos 
produced 9637 beam events. The processing of these events, particularly the scanning of 
emulsion films, is continuously going on. The 2011 run started on May 2011 and is still 
in progress. 

At the CNGS energies the average t~ decay length is submillimctric, so OPERA uses 
nuclear emulsion films as high precision tracking device in order to be able to detect such 
short decays. Emulsion films are interspaced with 1mm thick lead plates, which act as 
neutrino target and form the largest part of the detector mass. This technique is called 
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Emulsion Cloud Chamber (ECC). It was successfully used to establish the first evidence 
for charm in cosmic rays interactions [11] and in the DONUT experiment [12] for the 
first direct observation of the Vt- To date, nine Vt CC interactions have been observed 
by DONUT produced by a fixed target 800 GeV proton beam configuration. 

3. The OPERA Detector 

OPERA is a hybrid detector made of two identical Super Modules (SMI and SM2), 
each one formed by a target section and a muon spectrometer as shown in Figure 1. 
Each target section is organized in 31 vertical "waZ^s", transverse to the beam direction. 
Walls are filled with "ECC bricks" with an overall mass of 1.25 kton. They arc followed 
by double layers of scintillator planes acting as Target Trackers (TT) that are used to 
locate neutrino interactions occurred within the target. A target brick consists of 56 lead 
plates of 1 mm thickness interleaved with 57 emulsion films. The lead plates serve as 
neutrino interaction target and the emulsion films as 3-dimensional tracking detectors 
providing track coordinates with a sub-micron accuracy and track angles with a few 
mrad accuracy. The material of a brick along the beam direction corresponds to about 
10 radiation length and 0.33 interaction length. The brick size is 10cm x 12.5cm x 8cm 
and its weight is about 8.3kg. 



SM1 SIVI2 




VETO Target area Magnet and PT PT+XPC BMS 

RPCs 

Fig. 1. - View of the OPERA detector; the neutrino beam enters from the left. Arrows show 
the position of detector components, the VETO planes, the target and TT, the drift tubes (PT) 
laid out along the XPC, the magnets and the RPC installed between the magnet iron slabs. 
The Brick Manipulator System (BMS) is partly shown. 

In order to reduce the emulsion scanning load. Changeable Sheets (CS) [13] film 
interfaces have been used. They consist in tightly packed doublets of emulsion films glued 
to the downstream face of each brick. Charged particles from a neutrino interaction in a 
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brick cross the CS and produce signals in the TT that allow the corresponding brick to 
be identified and extracted by an automated Brick Manipulator System (BMS). 

The spectrometers consist of a dipolar magnet instrumented with active detectors, 
planes of RPCs (Internal Tracker, IT) and drift tubes (Precision Tracker, FT). Tasks of 
the spectrometers are muon identification and charge measurement in order to minimize 
the background. For muon momenta between 2.5 GeV/c and 45 GeV/c, the fraction of 
events with wrong charge determination is 1.2%. The to events ratio, within the 
selected momentum range, obtained from data can be directly compared with predictions 
based on Monte Carlo simulations: 3.92 ± 0.37(stoi.)% for data, 3.63 ± Q.\i{stat.)% for 
MC. Figure 2 left-side shows the momentum and momentum times charge distribution 
for data and MC. 
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p X charge (GeV/c) Reconstructed Bjorken-y 



Fig. 2. - Right: Muon charge comparison (momentum x charge): data (black dots with error 
bars) and MC (solid line) are normalised to one. Left: Bjorken-y variable reconstructed in data 
(dots with error bars) and MC (shaded areas). The MC distributions are normalised to data. 
The different contributions of the MC are shown in different colours: QE -I- RES contribution 
in light grey, DIS contribution in grey and the NC contamination in dark grey. 

In Figure 2 right-side, Bjorken-y distribution is shown for the events with at least a 
muon track. The agreement between data and MC simulation is reasonable. The sum 
of the QE and RES processes can be clearly seen as a peak at low y values. The NC 
contribution shows up at values of Bjorken-y close to one. The NC contribution becomes 
negligible when a track with its momentum measured by the spectrometer is required. 

A detailed description of the complete detector can be found in [7] . Event reconstruc- 
tion procedures and a performances of the OPERA electronic detectors can be found in 
more detail in [14]. 

4. Neutrino interaction location 

Neutrino event analysis starts with the pattern recognition in the electronic detectors. 
Charged particle tracks produced in a neutrino interaction generate signals in the TT 
and in the muon spectrometer. A brick finding algorithm is applied in order to select the 
brick which has the maximum probability to contain the neutrino interaction. The brick 
with the highest probability is extracted from the detector for analysis. The efficiency of 
this procedure reaches 83% in a subsample where up to 4 bricks per event were processed. 

After extraction of the brick predicted by the electronic detectors, its validation comes 
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from the analysis of the CS films. The measurement of emulsion films is performed 
through high-speed automated microscopes [15, 16] with a sub-micrometric position res- 
olution and angular resolution of the order of one milliradian. If no expected charged 
track related to the event is found in the CS, the brick is returned back to the detector 
with another CS doublet attached. If any track originating from the interaction is de- 
tected in the CS, the brick is exposed to cosmic rays (for alignment purposes) and then 
depacked. The emulsion films are developed and sent to the scanning laboratories of the 
Collaboration for event location studies and decay search analysis. 

All the track information of the CS is then used for a precise prediction of the tracks in 
the most downstream films of the brick (with an accuracy of about lOO/xm). When found 
in this films, tracks are followed upstream from film to film. The scan-back procedure is 
stopped when no track candidate is found in three consecutive films and the lead plate 
just upstream the last detected track segment is defined as the vertex plate. In order 
to study the located vertices and reconstruct the events, a general scanning volume is 
defined with a transverse area of 1 x \cm? for 5 films upstream and 10 films downstream 
of the stopping point. All track segments in this volume are collected and analysed. 
After rejection of the passing through tracks related to cosmic rays and of the tracks due 
to low energy particles, the tracks produced by the neutrino interaction can be selected 
and reconstructed. 

The present overall location efficiency averaged over NC and CC events, from the 
electronic detector predictions down to the vertex confirmation, is about 60%. 

5. Decay Search 

Once the neutrino interaction is located, a decay search procedure is applied to detect 
possible decay or interaction topologies on tracks attached to the primary vertex. The 
main signature of a secondary vertex (decay or nuclear inetaraction) is the observation 
of a track with a significant impact parameter (IP) relative to the neutrino interaction 
vertex. The IP of primary tracks is smaller than 10/im after excluding tracks produced by 
low momentum particles. When secondary vertices arc found in the event, a kinematical 
analysis is performed, using particle angles and momenta measured in the emulsion 
films. For charged particles up to about 6 GeV/c, momenta can be determined using 
the angular deviations produced by Multiple Coulomb Scattering (MCS) of tracks in the 
lead plates [17] with a resolution better than 22%. For higher momentum particles, the 
measurement is based on the position deviations. The resolution is better than 33% on 
1/p up to 12 GeV/c for particles passing through an entire brick. 

A 7-ray search is performed in the whole scanned volume by checking all tracks having 
an IP with respect to the primary or secondary vertices lower than 800/^m. The angular 
acceptance is ±500 mrad. The 7-ray energy is estimated by a Neural Network algorithm 
that uses the number of segments, the shape of the electromagnetic shower and also the 
MCS of the leading tracks. 

6. Data analysis 

In the following, the analysis results [18] of about 35% of the 2008 and 2009 data 
sample, corresponding to the 1.89 x 10^^ p.o.t are presented. The decay search procedure 
was applied to a sample of 1088 events of which 901 were classified as CC interactions. 
In the sample of CC interactions, 20 charm decay candidates were observed, in good 
agreement with the expectations from the Monte Carlo simulation, 16 ± 2.9. Out of 
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them 3 have a 1-prong topology where 0.8 ± 0.2 was expected. The background for the 
total charm sample is about 2 events. Several z^e-induced events have also been observed. 

Moreover, a first CC Vr candidate has been detected. The expected number of v-r 
events detected in the analysed sample is about 0.54±0.13(syst.) at Am^23 ~ 2.5 x 10^'^ 
eV"^ and full mixing. 

7. The first tau neutrino candidate 

In this section, the first tau neutrino candidate [18] will be described. The location 
and decay search procedure yielded a neutrino interaction vertex with 7 tracks. One 
track exhibits a visible kink with an angular change of 41 ± 2mrad after a path length 
of 1335 ± 35^171. The kink daughter momentum is estimated to be 12~t^ GeV/c by MCS 
measurement and its transverse momentum to the parent direction is 470l^2o MeV/c. 
The event is displayed in Figures 3 and 4. 




Fig. 4. Longitudinal view of the i^t candidate event. 
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Table I. - Kinematical variables of Vt candidate event. 



Variable Measured Selection criteria 



Kink angle (mrad) 


42 ±2 


>20 


Decay length (/^m) 


1335 ± 35 


Within 2 plates 


P daughter (GeV/c) 




>2 


PT daughter (MeV/c) 




>300 (7 attached) 


Missing PT {MeV/c) 


570^??° 


<1000 


Angle (j> (deg) 


173 ± 2 


>90 



All the tracks from the neutrino interaction vertex were followed until they stop or 
interact. The probability that one of them is left by a muon is estimated to be less than 
10~^. The residual probability for being a ly^CC event, with a possibly undetected large 
angle ji track, is about 1%; a nominal value of 5% is assumed. None of the tracks is 
compatible with being an electron. 

Two electromagnetic showers caused by 7-rays, associated with the event, have been 
located and studied. The energy of 7I is {5.6±1.0{stat.)±1.7{syst.)) GeV and it is clearly 
pointing to the decay vertex. The 72 has an energy of 1.2 ± OA{stat.) ± OA{syst.) GeV 
and it is compatible with pointing to either vertex, with a significantly larger probability 
to the decay vertex. 

All the selection cuts used in the analysis were those described in detail in the ex- 
periment proposal [19] and its addendum [20]. All the kinematical variables of the event 
and the cut applied are given in Table L 

The invariant mass of the two observed 7-rays is 120 ± 20{stat) ± 35{syst) supporting 
the hypothesis that they are emitted in a 7r° decay. The invariant mass of the charged de- 
cay daughter assumed to be a tt" and of the two 7-rays amount to 640lgQ^(stat)ltgQ°(syst) 
MeV/c, which is compatible with the p(770) mass. So the decay mode of the candidate 
is consistent with the hypothesis — > p^i^r (where the branching ratio is about 25%). 

8. Background Estimation 

The two main sources of background to the t~ — > h{mT'^)vr channel where a similar 
final state may be produced are: 

• the decays of charmed particles produced in j/^ CC interactions where the primary 
muon is not identified as well as the cc pair production in NC interactions 
where one charm particle is not identified and the other decays to a 1-prong hadron 
channel; 

• the 1-prong inelastic interactions of primary hadrons produced in v^GG interactions 
where the primary muon is not identified or in i/^ NC interactions and in which no 
nuclear fragment can be associated with the secondary interaction. 

The Monte Carlo expectation of the first background source is 0.007 ±0.004(syst.) event, 
the fraction produced in CC interactions is less than 10"'^ events. The second type of 
background amounts to 0.011 ± 0.006(s?/si.) event. The total background in the decay 
channel to a single charged hadron is 0.018±0.007(s?/si) events. The probability that this 
background events fiuctuate to one event is 1.8% (2.36(t). As the search for t~ decays 
is extended to all four channels, the total background then becomes 0.045 ± 0.023(s?/s<). 
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The probability that this expected background to all searched decay channels of the t~ 
fluctuates to one event is 4.5% {2.01a). At Arn^ = 2.5 x 10~^ eV"^ and full mixing, the 
expected number of observed r~ events with the present analyzed statistics is 0.54 ± 
0.13{syst) of which 0.16 ± 0.04(sys<) in the one-prong hadron topology, compatible with 
the observation of one event. 

9. Conclusions 

During 2008, 2009 and 2010 runs, a total intensity of 1.78 x 10^^, 3.52 x 10^^ and 
4.04 X 10^^ p.o.t. respectively, was achieved. Within these three years, 9637 beam events 
have been collected within the OPERA target. The neutrino interaction location and 
decay search are going on. 

A first candidate i^t CC interaction in the OPERA detector at LNGS was detected 
after analysis of a sample of events corresponding to 1.89 x 10^^ p.o.t. in the CERN CNGS 

beam. The expected number of Vr events in the analysed sample is 0.54±0.13(syst.). 
The candidate event passes all selection criteria, it is assumed to be a r~ lepton decaying 
into h^{mT'^)vr- The observation of one possible tau candidate in the decay channel 
h^(Tr'^)h'r has a significance of 2.36(t of not being a background fluctuation. 
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